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MULTI- WAVELENGTH OBSERVATIONS OF THE RADIO MAGNETAR PSR J 1622-4950 AND DISCOVERY OF ITS 

POSSIBLY ASSOCIATED SUPERNOVA REMNANT 
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ABSTRACT 

We present multi-wavelength observations of the radio magnetar PSR J 1622-4950 and its environment. 
Observations of PSR J1622-4950 with Chandra (in 2007 and 2009) and XMM (in 2011) show that the X- 
ray flux of PSR J 1622^4950 has decreased by a factor of ~ 50 over 3.7 years, decaying exponentially with a 
characteristic time of r = 360 ±11 days. This behavior identifies PSR J1622^4950 as a possible addition to the 
small class of transient magnetars. The X-ray decay likely indicates that PSR J1622^4950 is recovering from 
an X-ray outburst that occurred earlier in 2007, before the 2007 Chandra observations. Observations with the 
Australia Telescope Compact Array show strong radio variability, including a possible radio flaring event at 
least one and a half years after the 2007 X-ray outburst that may be a direct result of this X-ray event. Radio 
observations with the Molonglo Observatory Synthesis Telescope reveal that PSR J 1622^4950 is 8' southeast of 
a diffuse radio arc, G333. 9+0.0, which appears non-thermal in nature and which could possibly be a previously 
undiscovered supernova remnant. If G333. 9+0.0 is a supernova remnant then the estimates of its size and age, 
combined with the close proximity and reasonable implied velocity of PSR J 1622^1-950, suggests that these 
two objects could be physically associated. 

Subject headings: ISM: individual(G333. 9+0.0) - pulsars: individual (PSR J1622-4950) - radio continuum: 
stars - stars: neutron - supernova remnants - X-rays: stars 
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1. INTRODUCTION 

"Magnetar" has become the commonly used term to describe 
the emerging class of rare, young and highly magnetized neu- 
tron stars (B > 10 14 G) referred to as Anomalous X-ray Pul- 
sars (AXPs) and Soft Gamma Repeaters (SGRs), with perhaps 
as many as 23 now detected. 21 Magnetars are primarily bright 
X-ray emitters, for which most of this high energy radiation is 
thought to be generated through the decay of their strong mag - 
netic fields (for a review on magnetars see Mereghetti 2008). 
While magnetars may be X-ray luminous, only three such 
sources have been detected at radio wavelengths: XTE J 1 8 1 0— 
197 (Halpern et al. 2005; Camilo et al. 2006), IE 1547.0- 
5408 (Camilo et al. 2007b) and most recently PSR J1622- 
4950 (Levin et al. 2010). 

PSR J 1622^-950 is unusual as it is the first magnetar to have 
been discovered by its pulsed radio emission (Levin et al. 
2010). PSR J1622^950 was detected as a 4.3s period ra- 
dio pulsar in the High Time Resolution Universe survey per- 
formed with the Parkes 64m telescope (Keith et al. 2010), and 
was then recovered in other archival radio data-sets (Levin 
et al. 2010). Levin et al. (2010) showed this pulsar to be 
very different from ordinary radio pulsars in that it displays 
large variations in flux over time, with time spans of inactiv- 
ity during which it is undetected. Its variable pulse profile and 
high inferred magnetic field strength show that it has similar 
properties to the only other two known radio magnetars, XTE 
J1810-197 and IE 1547.0-5408. Both of these radio mag- 
netars also fit into an emerging subgroup known as transient 
magnetars, which randomly undergo bright X-ray bursts with 
increases in X-ray flux up to a factor of several 100 (for exam- 
ple see Tam et al. 2006; Bernardini et al. 2009, 201 1; Scholz 
& Kaspi 201 1). It appears that for both XTE J1810-197 and 
IE 1547.0-5408, the pulsed radio emission turned on as a re- 
sult of an X-ray outburst. 

In this paper we present multi-wavelength data that further 
confirm the magnetar identification of PSR J 1622^-950 and 
that support its transient nature. New and archival radio, 
infrared and X-ray observations are presented in Sections 2 
and the corresponding results are presented in Section 3. In 
Section 4 we use these data to explore PSR J1622-4950's 
X-ray and radio variability over the last ten years, and to 
study the variable nature of its polarized radio emission. We 
also discuss a possible 7-ray counterpart to PSR J 1622-4950, 
and identify a possibly associated young supernova remnant 
(SNR). 

2. OBSERVATIONS 
2.1. X-ray Observations 

2.1.1. Chandra Observations 

The position of PSR J 1622^-950 was observed three times 
with the Chandra X-ray Observatory, using the Advanced 
CCD Imaging Spectrometer (ACIS; Garmire et al. 2003). 
The first observation took place on 2007 June 13 as part of 
the "ChlcAGO" (Chasing the Identification of ASCA Galactic 
Objects) project. ChlcAGO is a survey designed to localize 
and classify the unidentified X-ray sources discovered during 
the ASCA Galactic Plane Survey (AGPS) (see Anderson et al. 

21 http://www.physics.mcgill.ca/~pulsar/magnetar/main.html 



2011, for some of the first ChlcAGO results; further details 
will be published by G. E. Anderson et al., in preparation.) 
PSR J1622^-950 was observed again on 2009 June 14 and on 
2009 July 10. These data were reduced using the Chandra In- 
teractive Analysis of Observations (CIAO) software, version 
4.3, following the online CIAO 4.3 Science Threads 22 . For 
details on these Chandra observations see Table 1 . 

2.1.2. XMM-Newton Observations 

We observed PSR J 1622^1950 with XMM-Newton starting 
on 2011 Feb 22. The PN (Struder et al. 2001) and MOS 
(Turner et al. 2001) cameras, which comprise the European 
Photon Imaging Camera (EPIC), were operated in Full Frame 
mode. Data were processed using SAS version 10.0.0, 23 and 
we employed the latest calibration files available at the time 
the reduction was performed (2011 March). Standard data- 
screening criteria were applied in the extraction of scientific 
products. The observation was highly affected by proton 
flares, which we had to cut from our data before proceed- 
ing with the scientific analysis, resulting in a net livetime 
exposure of 46.4 ks. Since both MOS and PN give consis- 
tent results, in the following we report only on the high-time- 
resolution PN data. Further details on the XMM EPIC-PN 
observation can be found in Table 1 . 

2.2. Infrared and Optical Observations 

We compared the Chandra position of PSR J1622-4950 (see 
Section 3.1.1. below) to optical and infrared point source 
catalogs and survey images including the US Naval Obser- 
vatory (USNO) B catalog, version 1.0 (Monet et al. 2003), 
the Two Micron All Sky Survey (2MASS; Skrutskie et al. 
2006), the Galactic Legacy Infrared Mid-Plane Survey Ex- 
traordinaire (GLIMPSE; Benjamin et al. 2003), and the 24 
/zm images from the MIPSGAL Survey (Carey et al. 2009) 
but no counterpart to PSR J 1622^-950 was identified in any 
of these data sets. On 2007 June 24 we observed the position 
of PSR J1622^1950 in the K s band for 13.5 minutes using 
the Persson's Auxiliary Nasmyth Infrared Camera (PANIC: 
Martini et al. 2004; Osip et al. 2008) on the 6.5m Magellan 
I (Baade) telescope, located at Las Campanas Observatory. 
The final image was a combination of 30s exposures carried 
out three times at each position of a nine-point dither pattern 
to account for the high background. Standard reduction pro- 
cedures were applied using the IRAF software package (Tody 
1986, 1993). We calibrated the field's photometry and astrom- 
etry using the source extraction software SExtractor (Bertin & 
Arnouts 1996) with comparisons to the Two Micron All Sky 
Survey point source catalog (2MASS PSC), which has a po- 
sition uncertainty of 0."1 (Skrutskie et al. 2006). No infrared 
counterpart was detected within 0."8 of the X-ray position of 
PSR J1622-4950 to a lower limit K s > 20.7. 

2.3. Radio Observations 
2.3.1. Archival 

PSR J 1622-4950 is coincident with a bright knot of emission 
that is part of a larger diffuse radio source seen in the first and 

22 http://cxc.harvard.edu/ciao/threads/ 

23 http://xmm.esac.esa.int/sas/ 
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second epoch Molonglo Galactic Plane Surveys (MGPS 1 and 
MGPS2 respectively; Green et al. 1999; Murphy et al. 2007) 
and also seen, albeit less clearly, in the continuum maps from 
the Southern Galactic Plane Survey (SGPS: Haverkorn et al. 
2006). This feature will be discussed in Section 3.2.1. The 
MGPS surveys were performed at 843 MHz with the Molon- 
glo Observatory Synthesis Telescope (MOST) at a resolution 
of 43" while the continuum SGPS images were created from 
Australia Telescope Compact Array (ATCA) 1 .4 GHz obser- 
vations at a resolution of 100". 

2.3.2. ATCA Observations 

We observed PSR J 1622^1950 with the ATCA simultane- 
ously at both 5 and 9 GHz, with a 128 MHz bandwidth in 
each band, on 2008 Nov 22 and 2008 Dec 5 as summarized in 
Table 2. A mosaic pattern was used in an attempt to both de- 
tect the magnetar and study the morphology of the surround- 
ing diffuse emission seen in the MGPS1, MGPS2 and SGPS. 
Both frequencies were observed for 6.79 and 8.01 hours in 
the EW367 and 750B configurations, respectively, using PKS 
B1934-638 for flux calibration and PMN J1603^1904 for 
phase and polarization calibration. Data reduction and anal- 
ysis were performed using the MIRIAD 24 software package 
using standard techniques. 

3. RESULTS 

3.1. X-ray Results 
3.1.1. X-ray Counterpart to PSR J1622^f950 

The 2007 June Chandra observation was an attempt to lo- 
calize the unidentified AGPS source AX J 162246^-946 as 
part of the ChlcAGO survey. Six point sources were de- 
tected within 4' of the position of AX J 162246^-946, all of 
which may have contributed to the X-ray emission originally 
detected with ASCA on 1999 April 22 as part of the AGPS. 
The brightest source, CXOU J162244.8-495054 (Evans et al. 
2010), which Levin et al. (2010) identified as the X-ray coun- 
terpart to PSR J 1 622-4950, was at least 20 times brighter than 
the other five sources. As CXOU J 162244. 8-495054 fell in 
the gap between CCDs for part of the 2007 June observation 
we only included those time intervals where the source was 
detected, resulting in an effective exposure time of ~ 2.02 ks 
and a count-rate of 96.8 ± 6.9 counts ks in the 0.3-8.0keV 
energy range. 

The 2009 July Chandra observation, two years later, also de- 
tected the same six sources in the 4' region surrounding the 
position of AX J162246-4946. While the 2009 July obser- 
vation was ^10 times longer than the 2007 June observa- 
tion, CXOU J162244.8^195054 had faded significantly with 
a much lower count-rate of 7. 1 ± 0.6 counts ks" 1 (equivalent 
to ~ 9.3 counts ks -1 when compared to the 2007 June ob- 
servation in the ACIS-S configuration 25 ) in the 0.3-8.0 keV 
energy range. The X-ray flux of CXOU J 162244.8^195054 
had therefore reduced by a factor of ~ 10 in the two years be- 
tween the 2007 June and 2009 July Chandra observations. 

24 http://www.atnf.csiro.au/computing/software/miriad/ 

25 All count-rate conversions are based on X-ray spectral fit of the 2007 
June Chandra observation, derived in Section 3.1.2., in order to account for 
the differences in the off-axis position and between the X-ray telescope in- 
struments. 



The best X-ray position for CXOU J 162244. 8^195054, to 
within a 0".8 radius circle for a 95% confidence, 26 is (J2000) 
RA = 16:22:44.89 and Dec = -49:50:52.7 taken from the 2009 
July Chandra observation. This agrees within 95% confi- 
dence of the X-ray position given by Levin etal. (2010). From 
now on we will refer to both the X-ray counterpart, CXOU 
J 162244. 8-495054, and radio counterpart of this magnetar as 
PSR J 1622^1950. 

The Chandra observation from approximately one month ear- 
lier (2009 June) also detected PSR J 1622^1950 very far off- 
axis, showing a count-rate slightly higher than that of the 2009 
July observation. This is equivalent to ~ 10.8 counts ks" 1 in 
the 0.3-8.0 keV energy range when compared to the 2007 
June Chandra observation. In 2011 Feb, our XMM observa- 
tion showed that PSR J 1622-4950 had faded even more, re- 
sulting in an EPIC-PN count-rate of 5.8 ±0.6 count ks" 1 in the 
0.4- 10.0 keV energy range (equivalent to ~ 1.8 counts ks" 1 
in the ACIS-S configuration of 2007 June in the 0.3 - 8.0 keV 
energy range). A summary of each X-ray observation can be 
found in Table 1 . 

It is also worth exploring to what extent the X-ray emis- 
sion from PSR J 1622-4950 may have contributed to the flux 
from AX J 162246^946 in the 1999 ASCA observation. PSR 
J 1622^1950 lies 4' away from the position of AX J 162246- 
4946 (Sugizaki et al. 2001). The other 5 point sources that 
were detected in both the 2007 June and 2009 July Chandra 
observations all fall within 3' of the position of AX J162246- 
4946. By fitting an absorbed power law to the spectrum of 
each of these point sources, assuming a power law spectral 
index of 2 and Nh = 2 x 10 22 cm" 2 , we estimate they con- 
tributed a total absorbed flux of ~ 1 x 10~ 13 erg cm" 2 s _1 in 
the 0.3-10.0 keV energy range. (The fluxes of 3 out of the 
5 point source were unchanged within a factor of 2 between 
the 2007 June and 2009 July Chandra observations. The flux 
of the fourth point sources may have decreased by a factor 
of ~ 4 between 2007 and 2009 and the flux of the fifth point 
source may have increased by a factor of ~ 6.) Using the 
2.0- 10.0 keV ASCA count-rates (Sugizaki et al. 2001) and 
the absorbed power law spectral fit described above, we esti- 
mate an absorbed X-ray flux of ~ 4 x 10~ 13 erg cm" 2 s" 1 from 
AX J 162246-4946. Subtracting the contribution from the 5 
nearby point sources described above, we estimate that the 
unabsorbed X-ray flux from PSR J 1622-4950 contributed, at 
most, ~ 3 x 10~ 13 erg cm" 2 s" 1 (corresponding to 75%) to the 
observed X-ray emission seen from AX J 162246-4946 in the 
AGPS. 

3. 1 .2. X-ray Spectrum and Variability 

In all three Chandra observations we extracted the X- 
ray spectrum of PSR J 1622^1950 using the CIAO tool 
specextract. As the source fell close to the edge of the 
CCD in the 2007 June observation and was very far off-axis 
in the 2009 June observation, specextract was run in ex- 
tended source mode to better handle the creation of response 
files. The spectra from the 2007 June and 2009 July Chan- 
dra observations were binned before fitting to include at least 
10 counts per bin, as this is the minimum number of counts 

26 This error takes into account the Chandra absolute positional accuracy 
(see http://cxc.harvard.edu/cal/ASPECT/celmon/) and the position uncer- 
tainties associated with the CIAO source detection algorithm wavdetect 
(Hong et al. 2005). 
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required to result in a statistically significant fit. The small 
number of counts detected, combined with the large absorp- 
tion in the Galactic plane, resulted in very few photons being 
detected below 2 keV in these two Chandra observations. A 
larger number of counts were detected during the 2009 June 
Chandra observation (the source had faded but the observa- 
tion was long with the source off-axis) so we chose to bin 
its spectrum to include at least 40 counts per bin. The XMM 
EPIC-PN source events and spectrum were extracted within a 
circular region of 20" centered on the peak of the point spread 
function of the source, while the background was obtained 
from a source-free region of similar size. When generating 
the EPIC-PN spectrum we have included events with PAT- 
TERN < 4 (i.e. single and double events). This XMM spec- 
trum was binned before fitting, using at least 40 counts per bin 
to compensate for the high background and not oversampling 
the intrinsic spectral resolution by more than a factor of 3. 27 

The X-ray spectra of magnetars are commonly fit using a 
blackbody plus power-law model or a multiple blackbody 
model. We chose to simultaneously fit the four X-ray ob- 
servations of PSR J 1622-4950 with both a single absorbed 
blackbody model and a single absorbed power-law model as 
the low number of counts preclude the identification of multi- 
ple components. The spectra were fit using XSPEC (Dorman 
& Arnaud 2001) with the absorption, Nh, set as a free pa- 
rameter but constrained to have the same value at each epoch, 
while the other parameters were set to vary individually. It 
should be noted that the resulting spectrum and fit to the June 
2007 Chandra observation of PSR J 1622-4950 may not be 
entirely accurate given its location close to the edge of the 
CCD, which resulted in the source dithering off the chip at 
regular intervals. By removing those time intervals where the 
source was not on the chip we have reduced this error but 
there can still be problems associated with the response files 
created during the extraction process. 28 

A simultaneous spectral fit with either an absorbed blackbody 
or an absorbed power-law model are equally reasonable for 
describing the four X-ray spectra of PSR J 1622-4950, both 
resulting in x 2 re d = 0-7. The blackbody and power-law spectral 
fits resulted in absorptions of Nh = 5.4+J Jj x 10 22 cm" 2 and 

Nh = 10.5^2 f x 10 22 cm" 2 (errors at 90% confidence), respec- 
tively, using abundances from Lodders (2003) and the pho- 
toelectric scattering section from Balucinska-Church & Mc- 
Cammon (1992) and Yan et al. (1998). These absorptions 
exceed the Galactic Nh at the position of PSR J1622aA§- 
4950 predicted from Hi surveys (Dickey & Lockman 1990; 
Kalberla et al. 2005) by a factor of 2-4, as expected at high 
column densities and low Galactic latitudes (e.g. Arabadjis 
& Bregman 1999). The absorbed blackbody and absorbed 
power-law fit parameters for each spectrum are listed in Ta- 
ble 3. 

In order to investigate the long-term X-ray variability of 
PSR J 1622-4950 we computed the absorbed and unabsorbed 
fluxes and their uncertainties for the four X-ray spectra, using 
the XSPEC model cf lux, in the 0.3- 10.0 keV energy range. 
These fluxes are listed in Table 3. (As the absorbed X-ray 

27 http://xmm.esa.int/sas/10.0.0/doc/specgroup/nodel4.html 
2S For caveats associated with the analy- 

sis of sources near the edge of a CCD see 
http://cxc.harvard.edU/ciao/threads/specextract/index.html#caveat.acisedge 
and http://cxc. harvard. edu/cgi-gen/wonderdesk/wonderdesk.cgi?do=faq_view 
_record&faq= 1 & view_detail= 1 &f aq_id=66 



fluxes calculated from the blackbody and power-law spectral 
fits are the same at each epoch within the 90% confidence we 
will use the fluxes obtained from the blackbody fit in the rest 
of our analysis unless otherwise stated. This investigation was 
also limited to the absorbed fluxes as the errors on the unab- 
sorbed power law fluxes are very large and unconstraining.) 
The 2007 June Chandra observation is our brightest X-ray 
detection of PSR J 1 622-4950, making it 9 times brighter than 
the two mid 2009 Chandra observations and 47 times brighter 
than in the 201 1 Feb XMM observation. As described in Sec- 
tion 3 . 1 . 1 . , PSR J 1 622^1950 has an X-ray flux upper limit of 
F x < 3 X 10" 13 erg cm" 2 s" 1 (0.3- 10.0 keV) in the 1999 April 
ASCA observation. This flux upper limit is > 6 times fainter 
than our X-ray flux measurement from the 2007 June Chan- 
dra observation but is consistent with the fluxes measured at 
the other three X-ray epochs. 

The best fit blackbody temperatures (kT) and power law spec- 
tral indices (T) of the simultaneous spectral fits are also listed 
in Table 3. While both of these parameters at each epoch 
are on the high end when compared to those seen from other 
known magnetars 29 they are not unreasonable (for example 
see Gelfand & Gaensler 2007; Enoto et al. 2009). The param- 
eters kT and V also indicate there was no significant spectral 
variability as these values remain the same between epochs 
within the 90% confidence. The spectra with the absorbed 
blackbodies are plotted in Figure 1 . 

We further investigated evidence of spectral evolution by 
comparing the hardness ratios resulting from simultaneously 
fitting all four spectra with an absorbed blackbody model and 
calculating the flux and 90% flux error for different energy 
bands (e.g. 0.3-2 keV, 2-4 keV, etc). The flux hardness ratios 
for different energy band combinations (soft/hard), propagat- 
ing through the 90% errors, shows no evidence of spectral 
variation between the epochs. While this method is model de- 
pendent it also demonstrates there is no strong evolution to 
PSR J1622-4950's X-ray spectrum. 

Levin et al. (2010) calculated a dispersion measure (DM) dis- 
tance of 9 kpc to PSR J 1622-4950 using the Cordes & Lazio 
(2002) NE2001 Galactic free electron density model. 30 We 
used this distance to estimate the unabsorbed 0.3- 10.0 keV 
luminosity at the four epochs from the unabsorbed fluxes cal- 
culated from both the blackbody and power-law spectral fits. 
These values are summarized in Table 3. The X-ray lumi- 
nosities calculated from the blackbody fit show a monotonic 
reduction over the four epochs and are consistent with the 
range of luminosities seen from other magnetars. 31 However, 
the luminosities calculated from the power-law fit are on the 
high end of known magnetar luminosities, with the large er- 
rors preventing us from observing any obvious evolution be- 
tween epochs. 

3.1.3. X-ray Timing Analysis 

Only the XMM data were used in our X-ray timing investiga- 
tion as our Chandra ACIS observations do not have sufficient 
time resolution to detect PSR J1622-4950's pulse period of 

29 See http://www.physics.mcgill.ca/~pulsar/magnetar/main.html 

30 This distance estimate should be treated with caution as proper motion 
studies have shown that DM-based distance measurements can differ by more 
than a factor of two to the distance calculated from parallax (Chatterjee et al. 
2009; Deller et al. 2009). 

31 See http://www.physics.mcgill.ca/~pulsar/magnetar/main.html 
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4.3s. Using the EPIC-PN observations, which have a time res- 
olution of 73.4 ms, we shifted all photons arrival times to the 
solar system barycenter, and used the XRONOS software for 
the timing analysis. We used data in the 0.3-4.0 keV energy 
range where the source was relatively bright, then folded the 
data using the pulsar ephemeris derived from ongoing radio 
monitoring: P= 4.32645312 s at Modified Julian Date (MJD) 
55586.5 (Levin et al., in preparation). We do not detect a sig- 
nificant signal in the X-ray band with a 3er upper limit on the 
0.3-4.0 keV pulsed fraction of 70% (derived as by Vaughan 
et al. 1994; Israel & Stella 1996). The pulsed fraction is de- 
fined as (N max -N min )/(N max +N min ), with N max and N min being 
the maximum and minimum counts of a putative sinusoidal 
signal at the frequency of PSR J 1622-4950. A similar anal- 
ysis was also conducted in the 0.3-2.0 and 2.0- 10.0 keV 
energy bands as the strength and phase of a magnetar's pul- 
sations can be energy dependent (Gotthelf & Halpern 2009). 
Each band had < 100 counts with no pulsations detected, re- 
sulting in an unconstraining pulse fraction upper limit of 99%. 

In order to search for short term variability from PSR J 1622- 
4950 in the three Chandra observations we corrected the pho- 
ton arrival times to the solar system barycenter. In the case 
of the 2007 June observation we also filtered out those times 
when the source was not on the chip to remove any contribu- 
tion caused by the source dithering on and off the chip. There 
was no evidence of short term variability beyond 40% of the 
mean count-rate between 250 and 2020s, 5000s and 60870s 
and 2500s and 20130s in the 2007 June, 2009 June and 2009 
July Chandra observations, respectively. 32 

3.2. Radio Results 

3.2.1. Archival 

The region surrounding PSR J 1622-4950 at radio wave- 
lengths in MGPS1 is shown in Figure 2(a). The position of 
PSR J 1622-4950 is denoted by a white "+" sign. The dif- 
fuse radio emission surrounding PSR J 1622-4950, encom- 
passed by the contours seen in Figures 2(b) and 2(c), is de- 
noted source A. There is a nearby Hll region, G333.6-0.2, 
located to the southwest of the pulsar (Paladini et al. 2003). 
The other three diffuse radio sources in the immediate vicin- 
ity of the pulsar are denoted B, C and D. Figures 2(b) and 2(c) 
show the corresponding infrared GLIMPSE 8 /im and MIPS- 
GAL 24 fim views of this area of sky, respectively, overlaid 
with MGPS 1 as black contours. 

The radio source A does not appear to have been catalogued 
in any surveys of this region. The brightest knot of radio emis- 
sion in source A is coincident with PSR J 1622^-950 and ex- 
tends out in some directions as far as 4'. This diffuse radio 
emission also forms a ring morphology below PSR J 1622- 
4950, which can be seen in Figure 2(a). This same mor- 
phology was also resolved by Levin et al. (2010) using the 
ATCA. Faint diffuse infrared emission is detected all over the 
immediate field as seen in Figures 2(b) and 2(c), but is not 
concentrated in the vicinity of the bright knot of radio emis- 
sion in source A. It is therefore more likely that this diffuse 
infrared emission originates from G333.6-0.2, which is ex- 
tremely bright in the radio and at 8 and 24 /im. At the South- 
ern edge of source A there is a very bright 24 /im source, 

32 The smaller number in each range is the size of a time bin chosen to 
include a statistically significant number of counts. The upper number is the 
exposure time. 



IRAS 16190-4946 (ellipse major and minor axis uncertainty 
of 23" and 3", respectively; Helou & Walker 1988), indicated 
by a white 23" radius circle in Figure 2(c). IRAS 16190-4946 
lies ~ 2' from PSR J 1622^-950 so is unlikely to be associated 
with source A immediately coincident with the magnetar. 

To the northwest of the pulsar there is a radio source in the 
shape of an arc, labeled B in Figure 2. Figures 2(b) and 
2(c) show no evidence for diffuse infrared emission match- 
ing the structure of source B. If we assume that this arc 
is part of a circular shell, we measure an equivalent radius 
of 7' ± 1' with approximate central coordinates of (J2000) 
RA=16:22:38 and Dec=-49:49:48. The inferred center of the 
arc is offset 0.8'-2.5' from the position of PSR J1622^I950. 
Another pulsar, PSR J 1622-4944, indicated by a black "+" 
sign in Figure 2, lies 1' inside the inner edge of the arc. The 
distance to PSR J1622^-944 is — 7.9 kpc calculated using 
the DM published by Manchester et al. (2001) assuming the 
Cordes & Lazio (2002) model. There is also an infrared dark 
cloud, SDC G333.900+0.022 (Peretto & Fuller 2009) coinci- 
dent with the arc but it is unlikely to be associated given that 
radio emission is not commonly associated with such clouds 
and that its morphology does not match that of source B. 

A possible counter-arc to source B, labeled source C in Fig- 
ure 2, is seen — 8' to the southeast of PSR J1622^-950 but 
the lack of symmetry between these arcs renders an associa- 
tion between the two speculative. Radio source D is coinci- 
dent with the 6.7 GHz maser source G333.93-0. 13 (Pestalozzi 
et al. 2005) and is likely also associated with IRAS 16194- 
4941 . At 24 /im other infrared sources are detected within the 
brighter radio contours of source D, which are very faint or 
undetected at 8 /zm. 

3.2.2. ATCA Results 

Initial analysis of the 2008 Nov and 2008 Dec ATCA data 
demonstrated that we did not have enough sensitivity to see 
the diffuse emission surrounding PSR J 1622-4950 as seen in 
MGPS, SGPS and the Levin et al. (2010) ATCA observations 
(which took place on 2009 December 8 and 2010 February 
27). While the Levin et al. (2010) ATCA observations used 
the same frequencies, array configurations, and integration 
times as the 2008 ATCA observations presented here, they 
were using the new ATCA correlator, which has a 2 GHz 
bandwidth and thus a much higher sensitivity. Regardless of 
the lack of diffuse emission, our ATCA observations did de- 
tect a radio point source at the magnetar's Chandra position. 
We analyzed the 2008 Nov and 2008 Dec observations sep- 
arately, using just the antenna 6 baselines (i.e., baselines be- 
tween 5 and 6 km). This significantly reduced sidelobe con- 
tamination from G333.6-0.2. 

The resulting time and phase averaged flux densities and po- 
larizations of PSR J 1622-4950, measured during the 2008 
Nov and 2008 Dec ATCA observations, are listed in Table 2. 
While the 9 GHz flux density remained fairly steady between 
the two observations, it increased by 22% at 5 GHz. The 
ATCA observations of PSR J 1622-4950 taken by Levin et al. 
(2010) (see Table 2) indicate a decrease in radio flux of ~ 68% 
and ~ 55%, at 5 and 9 GHz respectively, approximately one 
year after our 2008 Dec ATCA observation. 

The flux changes we observed between our 2008 Nov and 
2008 Dec ATCA observations demonstrated a steepening in 
the radio spectral index of PSR J 1622^-950 over two weeks, 
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from a = -0.13 ± 0.04 on 2008 Nov 22 to a = -0.44 ± 0.04 
on 2008 Dec 5, where S„ oc v a . These values contrast to the 
positive time-averaged spectral index derived by Levin et al. 
(2010), calculated using a combination of observations with 
different telescopes in the frequency range 1 .4 - 9 GHz and 
over various epochs between 1998 Feb and 2010 Jan. 

The polarization of PSR J 1622^-950 also changed signif- 
icantly between our two ATCA epochs. During the 2008 
Nov ATCA observation the linearly polarized fraction was 
~ 80%, and circular polarization < 8% at both frequencies. 
Two weeks later the linearly polarized fraction had dropped 
to < 20% but the source had become ^15% circularly polar- 
ized. (This circular polarization was negatively handed using 
the pulsar astronomy sign convention.) It should be noted that 
these polarization fractions quoted here are lower limits on 
the true peak polarization of PSR J 1622-4950 given that the 
ATCA observations of the pulse signal are time and phased 
averaged. 

4. DISCUSSION 

4. 1 . Variability 

The spin-down luminosity of PSR J 1622^-950 calculated by 
Levin et al. (2010) is E = 8.5 x 10 33 erg s" 1 . Using the unab- 
sorbed blackbody luminosities listed in Table 3 we find that 
L x ~ 3.5E, L x ~ 0.6£, L x ~ 0.5£ and L x ~ 0.1E for the 2007 
June, 2009 June, 2009 July Chandra and 201 1 Feb XMM ob- 
servations, respectively, in the 0.3- 10.0 keV energy range. 
All these ratios are significantly higher than what we expect 
from young pulsars but are similar to what we see from other 
magnetars. One of the defining magnetar characteristics is 
L x > E (Mereghetti 2008) but this relationship can also be the 
case for cooling neutron stars (for example see van Kerkwijk 
& Kaplan 2008; Halpern & Gotthelf 201 1). However, the X- 
ray variability and relatively high X-ray luminosity that we 
observe from PSR J 1622^-950 are not consistent with cool- 
ing. We therefore confirm the Levin et al. (2010) identifica- 
tion of PSR J 1622-4950 as a magnetar. 

The simultaneous blackbody and power-law spectral fits to 
our four X-ray spectra of PSR J 1622^-950 show that over 
3.7 years the X-ray flux decreased by a factor of 47. Fig- 
ure 3 is a light-curve showing the blackbody absorbed X-ray 
flux values in blue where the Chandra observations are de- 
noted by open circles and the XMM observation is a filled 
circle. Figure 3 clearly illustrates the smooth decay of the X- 
ray emission from an outburst that may have occurred before 
or during the 2007 June Chandra observation. We therefore 
suggest that PSR J 1622^-950 is a transient magnetar, having 
demonstrated similar variations in X-ray flux to the two tran- 
sient magnetars XTE J1810-197 and IE 1547.0-5408 (Got- 
thelf & Halpern 2005, 2007; Bernardini et al. 2009; Halpern 
et al. 2008). 

The decays of outbursts from transient magnetars have been 
modeled as exponentials, power laws and multiple power laws 
(for example see Figure 3 of Rea & Esposito 201 1). We fitted 
both an exponential and power-law to the absorbed blackbody 
X-ray flux values from the Chandra and XMM observations. 
(We determined from these X-ray observations that the cross - 

33 Typical E and L x values for most magnetars are summarized at 
http://www.physics.mcgill.ca/~pulsar/magnetar/main.html 



calibration between the instruments was no larger than 12%, 
which is within the 90% confidence flux errors.) The decay 
of the X-ray light-curve is best described by an exponential 
with a characteristic decay time of r ~ 360 ±11 days. This is 
a similar to the X-ray decay of XTE J1810— 197 (Bernardini 
et al. 2009) but much slower than the decay of the multiple 
bursts produced by IE 1547.0-5408 (for example see Israel 
etal. 2010). 

In Section 3.1.2. we demonstrated that there is no evidence of 
strong evolution in PSR J1622^1950's X-ray spectrum over 
the 3.7 years of X-ray observations. A consistent temper- 
ature between epochs was also observed from XTE J1810- 
197, where the blackbody temperature remained unchanged 
for over three years after the initial outburst (Bernardini et al. 
2009). Further X-ray observations with higher count statis- 
tics are required to determine if the X-ray spectrum of PSR 
J 1622^-950 is evolving with time. 

If the fading X-ray emission that we have observed from 
PSR J 1622-4950 is thermal, then it is likely emitted by a 
hot spot on the surface of the magnetar that remains constant 
in temperature but decreases in size over time. This behav- 
ior is potentially explained by the magnetar coronal model 
(Beloborodov & Thompson 2007). Specifically Gotthelf & 
Halpern (2005) suggest that the shrinking of a hot spot after a 
burst could be due to the decay of currents or rearrangement 
of the magnetic field lines altering the heat being channeled to 
the surface. The changes in X-ray flux on the order of weeks 
to years may be the magnetar's crust plastically responding to 
the unwinding of these fields, which in turn deposits energy 
into the magnetosphere resulting in transient behavior (Muno 
et al. 2007, and references therein). 

The radio spectral index of PSR J 1622-4950 (listed in Ta- 
ble 2) appears to be variable and far flatter than the steep, 
a « -1.6, stable average spectral index we expect from or- 
dinary young radio pulsars (Lorimer et al. 1995). It is, how- 
ever, very similar in behavior to the other radio magnetars, 
XTE J1810-197 and IE 1547.0-5408 (Camilo et al. 2007c; 
Lazaridis et al. 2008; Camilo et al. 2008). The time-variable 
radio spectral indices of magnetars are not well understood 
but Thompson (2008a,b) suggests that this phenomenon could 
be the result of current-driven instabilities in the closed mag- 
netosphere of the magnetar. Serylak et al. (2009) (and refer- 
ences therein) also speculate that the observed flat radio spec- 
tra of magnetars could be the result of the open magnetic field 
lines having a high plasma density. 

As mentioned in Section 2.2. our PANIC observations ob- 
tained a counterpart lower limit of K s > 20.7 for PSR J1622- 
4950. Currently only seven magnetars have been detected in 
the ZQ-band (Rea & Esposito 201 1). By correcting for the dif- 
ference in extinction and relative distances, a /Q-band coun- 
terpart similar to that seen from the majority of these infrared 
magnetars, at the position of PSR J 1622-4950, would be ex- 
tremely faint, having a magnitude > 24. 7. 34 While the vari- 
able infrared behavior is not consistent between magnetars the 
ratio between their X-ray and infrared flux appears to remain 
fairly consistent where F v v ~ 1 x 10 4 /^ (Durant & van Kerk- 
wijk 2005), predicting a TT-band magnitude of 24.5 for PSR 



34 For details on the X'-band observations and distance 
estimates of each magnetar see Rea & Esposito (2011), 
http://www.physics.mcgill.ca/~pulsar/magnetar/main.html and references 
therein. 
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J 1622^-950 based on the 2007 June Chandra observation. A 
much deeper infrared observation is therefore required to de- 
tect the counterpart to PSR J1622-4950 during a future period 
of X-ray bursting activity. 



4.2. Polarization 

In Section 3.2.2. we showed that the linearly and circularly 
polarized fractions of the radio emission from PSR J 1622- 
4950 both changed significantly between 2008 Nov and 2008 
Dec. In contrast, phase resolved radio observations of the 
two radio magnetars XTE J1810-197 and IE 1547.0-5408 
show that their linear polarization is consistently high and that 
neither their linear or circular polarization exhibit dramatic 
changes in intensity over time (Kramer et al. 2007; Camilo 
et al. 2007d, 2008). 

The significant linear and circular polarization variability we 
observed from PSR J 1622-4950 can be explained by changes 
in the geometry of the magnetosphere of the magnetar caus- 
ing the overall pulse profile to vary (Camilo et al. 2007a). 
The 2008 Nov and 2008 Dec ATCA observations were phase- 
and time-averaged, amplifying the observed polarization vari- 
ability over the intrinsic behavior. Results from Levin et al. 
(2010) show that, unlike most normal radio pulsars, the ra- 
dio pulse profile of PSR J 1 622^-950 is variable between con- 
secutive epochs, similar to the radio pulse profiles of XTE 
J1810-197 and IE 1547.0-5408 (Camilo et al. 2007a; Kramer 
et al. 2007; Camilo et al. 2008). In Figure 1 of Levin et al. 
(2010), we see that in three consecutive epochs, the pulse 
profile shape changed from a double-peaked profile to one in 
which only the leading peak was detected. Figure 4 of Levin 
et al. (2010) shows that over the entire pulse profile at 1 .4 and 
3.1 GHz, when both peaks were detected, the polarized po- 
sition angle (PPA) swung through 180 degrees. If an ATCA 
observation of PSR J 1622^-950 takes place when the entire 
double-peaked pulse is switched on, and is therefore experi- 
encing a very large PPA swing, phase-averaging of the pulse 
would produce a low linear polarization measurement, such 
as what we observed in 2008 Dec. However, if only part of 
the pulse is switched on, such as in the case of the last pulse 
profile shown in Figure 1 of Levin et al. (2010), then a phased- 
averaged observation does not experience the full PPA swing. 
The phase-averaging would then result in less depolarization 
such as in the case of our 2008 Nov ATCA observation. 

The difference in circular polarization between the two ATCA 
epochs can be explained by the same changes in the pulse 
profile. In 2008 Dec, when we assume that the entire pulse 
profile is switched on, we see a significant fraction of cir- 
cular polarization that was not observed in 2008 Nov. It is 
possible that the part of the pulse profile that was switched 
off during the 2008 Nov observation was the circularly po- 
larized component, resulting in a lack of detectable circularly 
polarized emission during that epoch. This component then 
switched back on again in 2008 Dec, allowing us to detected 
circular polarization. Observations by Levin et al. (2010) have 
demonstrated that changes in the pulse profile of PSR J 1622- 
4950 happen on the timescale of days and could therefore 
have occurred in the 14 days between the 2008 Nov and 2008 
Dec ATCA observations, resulting in the observed polariza- 
tion variability. 



4.3. Correlation between Time Variability in Radio and 
X-rays 

In order to determine if there is any correlation between the X- 
ray and radio emission from PSR J 1622-4950, we have com- 
pared our X-ray and ATCA results with the Parkes 1 .4 GHz 
light-curve reported by Levin et al. (2010). Figure 3 shows 
the light-curve of PSR J 1622-4950 where the radio 1.4 GHz 
flux values are in red (flux scale on the left axis) and X-ray 
flux values are in blue (flux scale on the right axis). The ra- 
dio flux points include the Parkes measurements from Levin 
et al. (2010) and our ATCA detections extrapolated to 1.4 
GHz. This extrapolation assumes that the radio spectral index 
of PSR J 1622^1950, calculated from these ATCA observa- 
tions (see Table 2), describes the radio spectrum down to low 
frequencies, just as in the case of XTE J1810— 197 (Lazaridis 
et al. 2008). This results in estimated 1.4 GHz fluxes of ~ 39 
mJy during the 2008 Nov observation and ~ 69 mJy during 
the 2008 Dec observation, which are far brighter than any 
other radio detection of PSR J 1622-4950 at 1.4 GHz. If our 
values represent the fluxes of PSR J 1622-4950 at 1.4 GHz 
on these two dates then this magnetar was undergoing a radio 
flaring event during our ATCA observations and was there- 
fore in a high radio state. All other radio measurements, in- 
cluding those at earlier epochs than covered by the time range 
of Figure 3 (see Figure 1 of Levin et al. 2010), have revealed 
considerably lower flux values or only upper limits. 

The 2007 June (MJD 54264) Chandra observation detected 
PSR J 1622-4950 in a reasonably high X-ray flux state when 
the magnetar was likely recovering from a recent X-ray out- 
burst. The 2008 Nov and 2008 Dec ATCA observations then 
indicate a high radio flux from PSR J 1622-4950 1.45 years 
after this Chandra observation (assuming that the spectral in- 
dex from the ATCA observations can be extrapolated to 1 .4 
GHz). Studies ofXTE J1810-197 and lEJ1547.0-5408show 
that their radio emission was trigged by an X-ray outburst 
(Halpern et al. 2005; Camilo et al. 2007b) and, in the case 
of XTE J1810— 197, this radio emission peaked in intensity 
~ 3 years after its X-ray outburst (Camilo et al. 2006). There- 
fore, it is possible that the high radio state observed for PSR 
J 1622-4950 with the ATCA in 2008 was triggered by the X- 
ray outburst that occurred around the time of the 2007 June 
Chandra observation. 

The 2009 June and 2009 July Chandra observations show that 
the X-ray flux declined by a factor of 8 following the 2007 
June Chandra observation. The 2009 June and 2009 July 
Chandra observations occurred during the period between 
MJD 54900 and MJD 55250 for which the Parkes observa- 
tions showed the 1 .4 GHz flux to be extremely variable (see 
Figure 3). The XMM observation shows that by 201 1 Feb the 
X-ray flux had declined by an additional factor of 6. Further 
X-ray and radio observations are required to determine if PS R 
J1622-4950 has subsequently returned to a quiescent state. 

An X-ray flux upper limit of F x < 3 x 10~ 13 erg cm" 2 s" 1 
(0.3-10.0 keV) was also obtained for PSR J1622-4950 in 
the 1999 April (MJD 51291) AGPS observation (see Section 
3.1.1.). Figure 1 of Levin et al. (2010) indicates that Parkes 
did not detect PSR J 1622-4950 at 1.4 GHz before, during or 
shortly after this ASCA observation. It was not until 50 days 
after the ASCA observation that Parkes detected radio emis- 
sion. Unfortunately, without further X-ray flux history, it is 
not possible to determine if the X-ray detection is real and/or 



G. E. Anderson et al. 



related to the subsequent radio detection. 

4.4. "/-ray Counterpart 

PSR J1622^I950 falls within the 95% error circle of the 
Fermi Large Area Telescope (LAT) source OFGL J1622.4- 
4945 (Abdo et al. 2009). OFGL J1622.4-^945 was one of the 
Fermi-LAT bright 7-ray sources (Abdo et al. 2009), detected 
in the first three months of observations, with a statistical sig- 
nificance > 10(7. However, OFGL J 1622.4^945 is one of 10 
Fermi-LAT Bright Source List sources that were not detected 
in the Fermi-LAT First Source Catalog (1FGL), a more recent 
catalog in which each source detection is based on the average 
flux over an eleven month period with a statistical significance 
> 4a (Abdo et al. 2010a). This is not unexpected due to its lo- 
cation in the Galactic ridge, since Fermi sources in this region 
are far more difficult to detect and characterize in the 1FGL 
catalog analysis (Abdoet al. 2010a). 

PSR J1622-1950 lies - 5.8' from the centroid of OFGL 
J 1622.4-4945 so it is worth investigating if there is an associa- 
tion between these two objects. Fermi 7-ray pulsars have been 
found to have spin-down luminosities between £ = 3x 10 33 
and 5 x 10 38 erg s" 1 (Abdo et al. 2010c). While the spin- 
down luminosity of PSR J 1622^-950 is encompassed in this 
range, this magnetar is far more distant than the 7-ray pulsars 
listed in Abdo et al. (2010c), the majority of which lie within 
3kpc of Earth. The Fermi 7-ray pulsars have an E/d 2 be- 
tween 3 x 10 33 and 1 x 10 38 erg s" 1 kpc" 2 (Abdo et al. 2010c). 
The value of E/d 2 for PSR J1622^950 is 1 x 10 32 erg s" 1 
kpc" 2 , an order of magnitude smaller than the minimum value 
seen from the Fermi pulsars. An association between OFGL 
J 1 622.4-4945 and PSR J 1 622-4950 is therefore unlikely, par- 
ticularly as an investigation of much closer magnetars did not 
yield a 7-ray detection with Fermi-LAT (Abdo et al. 2010b). 

4.5. Supernova Remnant Association 

Magnetars are young neutron stars. However, while the 
number of identified magnetars is increasing, few have been 
found to have a convincing association with a supernova rem- 
nant (SNR): IE 2259+586 in CTB 109 (Fahlman & Gre- 
gory 1981), IE 1841-045 in Kes 73 (Vasisht & Gotthelf 
1997), IE 1547.0-5408 in SNR G327.24-0.13 (Gelfand & 
Gaensler 2007) and the yet to be confirmed magnetars AX 
J1845-0258 in SNR G29.6+0.1 (Gaensler et al. 1999) and 
CXOU J171405.7-381031 in CTB 37B (Aharonian et al. 
2008; Halpern & Gotthelf 2010). 35 Such associations are im- 
portant because SNRs provide independent constraints on the 
environment and the properties of the associated magnetar. 

Source B (and its possibly associated counter-arc source C) 
appears to have similar properties to SNRs as the apparent 
lack of associated diffuse infrared emission, as mentioned in 
Section 3.2.1., implies a non-thermal nature (for example see 
Brogan et al. 2006). PSR J1622^950 also resides within 2.5' 
of the center of the extrapolated structure. Another pulsar, 
PSR J 1622^1944, mentioned in Section 3.2.1., lies ~ 1' from 
Source B but has a high characteristic age of 10 6 yrs (Manch- 
ester et al. 2001). While it is true that a pulsar's characteristic 

35 Another possible association is SGR 0526-66 with SNR N49 in the 
Large Magellanic Cloud (Cline et al. 1982; Gaensler et al. 2001; Park et al. 
2012). 



age can be overestimated, it is extremely unusual for a pulsar 
with such a high characteristic age to be found associated with 
an SNR. In the absence of any further evidence, PSR J 1622- 
4944 appears to be an ordinary old radio pulsar unrelated to 
source B. 

Confirmation of source B as an SNR requires the measure- 
ment of a non-thermal radio spectral index or the detection 
of linear polarization, neither of which is possible to obtain 
with our current data-sets. (SGPS lacks the sensitivity re- 
quired to disentangle the diffuse sources in this complicated 
region and the MGPS does not have any polarization infor- 
mation. The frequencies of both these surveys are also too 
closely spaced to provide a meaningful spectral index esti- 
mation.) Instead we need to consider the feasibility of such 
an identification by exploring a possible connection between 
source B and PSR J 1622-4950 using some of the pulsar/SNR 
association criteria established by Kaspi (1996). In the fol- 
lowing discussion the criteria we explore are whether PSR 
J1622^I950 and G333.9+0.0 have consistent ages and if the 
implied transverse velocity of the magnetar away from the 
assumed explosion site is reasonable. As we do not know the 
distance to G333. 9+0.0 we cannot investigate whether the dis- 
tances to the magnetar and SNR are consistent. For the pur- 
poses of further consideration we assume a common distance 
of 9 kpc to both sources. The probability of chance alignment 
between the magnetar and SNR also needs to be considered as 
demonstrated by Gaensler et al. (2001). We designate source 
B as G333.9+0.0 based on its approximate centroid and will 
assume it is an SNR in our discussion below. 

The positional coincidence of PSR J 1622-4950 with the cen- 
ter of G333. 9+0.0 could suggest a possible association. The 
chance probability of finding an arc in MGPS 1 whose emis- 
sion is non-thermal (based on comparisons with GLIMPSE 
data), and whose center is within 2.5' of a given position on 
the sky, is about 5% based upon inspecting 100 random posi- 
tions for 315 < I < 357 and \b\ < 0.4. This probability is not 
particularly high or low, and so does not strengthen or argue 
against an association. 

The angular separation between the estimated center of the 
arc G333.9+0.0 and the position of the PSR J1622-4950 is 
0.8-2.5'. Using PSR J 1622^1950's characteristic age of 4 kyr 
(Levin et al. 2010), this leads to a projected velocity of 500- 
1500 km s" 1 for a distance of 9 kpc. This is consistent with the 
overall velocity distribution for pulsars (Arzoumanian et al. 
2002), but higher than that observed or inferred for magnetars 
(Gaensler et al. 2001; Helfand et al. 2007; Kaplan et al. 2009; 
Deller et al. 2012; Tendulkar et al. 2012). 

No X-ray emission is detected from G333.9+0.0 in the XMM 
observation; the upper limit on the count-rate is 0.04 counts 
s" 1 using standard errors and roughly accounting for vi- 
gnetting (Romer et al. 2001). Assuming a thin-shell mor- 
phology and the Sedov-Taylor solution (Sedov 1946a,b; Tay- 
lor 1950a,b), we find that the upper limit implies a preshock 
ambient density lower than «o = 0.05 cm" 3 for an explosion 
energy of 10 51 erg. For a SNR with a 18 pc radius this am- 
bient density predicts a swept up mass of 30 M Q , support- 
ing Sedov-Taylor expansion. Here we have used the stan- 
dard shock jump conditions for an ideal gas to determine the 
gas temperature based on the Sedov-Taylor age for the given 
value of no, and then estimated the X-ray count rate assuming 
a Raymond-Smith plasma model in XSPEC. We find a Sedov- 
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Taylor upper-limit age of 6 kyr, at most a factor of 1 .5 greater 
than PSR J1622^950's characteristic age of 4 kyr, making 
G333. 9+0.0 potentially young and similar in age to convinc- 
ing magnetar/SNR associations (Gaensler et al. 2001). This 
Sedov-Taylor age leads to an upper limit on shock velocity 
of ~ 1200 km s . The associated proton temperature is < 2 
keV, indicating an electron temperature of < 0.2 keV using 
the electron-ion equilibration relation from Ghavamian et al. 
(2007). While the above upper limit on the density is much 
lower than the mean ISM value, it is not unreasonable for the 
low-density cavities formed from stellar winds of the massive 
stars that lead to core-collapse supernovae. We conclude that 
the radius of the shell and the absence of X-ray emission are 
both consistent with young SNRs in a low density medium. 

The implied non-thermal nature of source A, combined with 
the positional coincidence of PSR J 1622-4950 within its ex- 
tent, raises the possibility that source A could be a pulsar wind 
nebula (PWN) generated by PSR J 1622^1950. This would 
be an interesting result, as there is currently no known radio 
PWN associated with a magnetar. However, a radial profile 
analysis, comparing the detection of PSR J 1622^-950 to the 
point-spread-function of the XMM image, do not show any 
evidence of an extended X-ray source. We also search for ex- 
tended X-ray emission resulting from a dust-scattering halo, 
similar to that seen around IE 1547.0-5408 (Tiengo et al. 
2010; Olausen et al. 2011), but none was detected. ATCA 
observations, with longer integration times, are required to 
determine the true nature of source A. 

The above do not provide direct evidence that G333. 9+0.0 is 
the shell of an SNR associated with PSR J 1622-4950, but 
there is no firm evidence to argue against such an association 
either. Another possible SNR association includes the ring of 
diffuse radio emission that forms part of source A, which sits 
- 2' south of PSR J 1622-4950. Levin et al. (2010) discuss 
the possibility that this ring could be the parent SNR to PSR 
J 1622-4950 as it appears non-thermal in nature given that it 
lacks an infrared counterpart. However, they consider a link 
unlikely given the ring's small size and the high implied mag- 
netar birth velocity. 

5. CONCLUSIONS 

In this paper we have confirmed the Levin et al. (2010) mag- 
netar identification of PSR J 1622^-950 through the detection 
of significant X-ray flux variability and high X-ray luminos- 
ity. The high dynamic range in the X-ray flux, combined 
with the exponential characteristic decay time of r = 360 ±11 
days, suggests that PSR J 1622^-950 may be a new addition 
to the transient magnetar class, and could possibly be recov- 
ering from an X-ray outburst that occurred before or during 
the 2007 June Chandra observation. This X-ray flux vari- 
ability, along with the variable radio flux and spectral index, 
make PSR J 1622-4950 similar to the two other known radio 
magnetars, XTE J1810-197 and IE 1547.0-5408. Observa- 
tions with the Australia Telescope Compact Array show that 
PSR J 1622-4950 may have undergone a radio flaring event 
approximately one and a half years after the 2007 June Chan- 
dra observation, which could have been triggered by the X- 
ray outburst that occurred around this time in 2007. The prox- 
imity of PSR J 1 622^-950 to the supernova remnant candidate 
G333. 9+0.0, the implied transverse velocities for PSR J1622- 
4950, and the apparent young age of SNR G333.9+0.0, all 
support the possibility of a new magnetar/SNR association. 



ASCA observations in 1999 indicate that PSR J 1622^1950 
may be responsible for some of the X-ray flux detected from 
AX J 162246^1956 in the ASCA Galactic Plane Survey (Sug- 
izaki et al. 2001). In the AGPS the magnetars XTE J1810- 
197, IE 1547.0-5408, IE 1841-045 (Vasisht & Gotthelf 
1997), and SGR 1806-20 (Laros et al. 1987; Ulmer et al. 
1993), were all detected as X-ray sources. Through our work 
in ChlcAGO we anticipate the discovery of other magnetars, 
which will allows us to further define the properties of this 
unusual population of neutron stars. 
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Table 1 

X-ray Observations of PSR J 1622^950 



Telescope 



Obs ID 



Date" 



MJD C 





Instrument 


Exp Time (ks) b 


Chandra 


8161 
ACIS-S 


2007-06-13 
2.02 


Chandra 


9911 
ACIS-I 


2009-06-14 
60.10 


Chandra 


10929 
ACIS-I 


2009-07-10 
19.90 


XMM 


0654110101 
EPIC-PN 


2011-02-22 
46.4 



Count-rate 



counts ks 



ACIS-S rate c 



counts ks 



Off-axis Angle 



54264 



54996 



55022 



55615 



96.8 ±6.9 

9.7 ±0.6 
7.1 ±0.6 

5.8 ±0.6 



10.8 

9.3 
1.8 



4.0 



16.1 



1.1 



a The date is in the form yyyy-mm-dd. 

The quoted exposure times are the effective exposure time after time intervals when there was flaring or when the source was on the chip gap have been removed. 
c Modified Julian Date 

Total observed count-rate in the 0.3 — 8.0 keV energy range. 
e Model predicted ACIS-S equivalent count-rate when compared to the 2007 June Chandra observation in the 0.3-8.0 keV energy range. 



Table 2 
ATCA Observations of PSR J 1622^4950 



Date 


MID 


Central Frequency 
(MHz) 


Flux Density 
(mly) 


Spectral Index 
(a) 




Polarization 




(yyyy-mm-dd) 


Linear (mJy) 


PA (deg) 


Circular (mJy) 


2008-11-22 


54793 


5312 
8768 


33.0 ± 0.3 
30.9 ± 0.6 


-0.13 ± 0.04 


26.6 ± 0.7 (79%) 
25.0 ±0.8 (81%) 


-17.5 ±0.5 
-25.8 ±0.7 


< 2.0 (< 6%) 

< 2.5 (< 8%) 


2008-12-05 


54806 


4800 
8256 


40.4 ± 0.3 
31.9 ±0.6 


-0.44 ± 0.04 


5.7 ±0.4 (14%) 

5. 8 ±0.7 (18%) 


+26.7 ±1.5 

-22.5 ±2.5 


-6.2± 0.3 (15%) 
-4.8± 0.5 (15%) 


2009-12-08 
2010-02-27 a 


55174 
55255 


5500 
9000 


13±1 
14.3 ±0.8 


+0.2 + 0.2 









Note. — All errors are 1 a. The circular polarization upper bounds for the 2008-1 1-22 observations are limits on the magnitude. 
' The flux densities quoted are the average values at 5.5 and 9 GHz from the 2009-12-08 and 2010-02-27 ATCA observations, taken from Levin et al. (2010). 



Table 3 

X-ray Spectral modeling of PSR J 1622^+950 



Telescope 




Absorbed Blackbody Fit" 






Absorbed Power-law Fit a 


Date 


kT 

Nh 


A- red 


Lx.unab 


i" 

Nh 


Fx.abs 
X-red 


A.v, it nab 

Lx.unab 


Chandra 
2007-06-13 


0.7 ±0.1 

5 4+1.6 


1.4 ± 0.3 x 10-' 2 
0.7 


3.5+> 4 x 10- 12 
3.4^; 4 X 10 34 


4 9+ 1 - 
10.5+|f 


1.5±°- 4 x 10-' 2 
0.7 


3.3 +24 7 4 x IO" 10 
■i.lf* x 10 36 


Chandra 
2009-06-14 


0.7 ±0.1 

5.4+^ 


1.8+°;3 X 10- 13 


4.9± 2 ' 2 x io- 13 
4.8+ 2 ^ x IO 33 


43+0.9 

la5 ll 


2.1+o ; 4 xlO-' 3 


5.9% l x ur 11 
5.7^ 8 8 x IO 35 


Chandra 
2009-07-10 


0.8$ 

J.+ , 4 


1.6! u: * x IO" 13 


3.6+JI x 10- 13 
3.5^;jj x IO 33 


3 6 +L2 
J -°-i.i 

10.5! 2 ;f 


2.0^1 x 10-' 3 


lJiPi 6 x 10" 11 

1.2+j 2 , 2 x IO 35 


XMM 

2011-02-22 


0.5 ±0.1 
5 4 +L6 


3-0+g | x 10- 14 


l.l+<" x 10- 13 

i.r«ft x io 33 


5 4 +L3 
-UK 


3.2^ x IO- 14 


7.9 + " 3 - 3 xl()-" 
l.T™ A x IO 35 



a The best fit absorbed blackbody and absorbed power-law model parameters including temperature, kT (keV), spectral index, V, absorption column density, Nh (10 cm ), and the 
reduced chi-square, X re d^ fr° m Chi Gehrels statistics. The absorbed and unabsorbed X-ray flux, F XM j ?s and F XMna b (erg cm" s~ ) respectively, as well as the unabsorbed X-ray luminosity, 
L x ,unab ( er g s )> are over the 0.3 - 10.0 keV energy range. The unabsorbed luminosity was calculated assuming a distance of 9 kpc (Levin et al. 2010). All fit parameter errors are for 

90% confidence. 
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Figure 1. The 2007 June (black), 2009 June (red) and 2009 July (green) Chandra spectra and 2011 Feb (blue) XMM spectrum of PSR J1622^1950. These 
spectra were fit simultaneously with an absorbed black body model with Nu locked between epochs but with kT and the normalization values for each spectra 
allowed to vary individually. The bottom panel shows the residuals of these fits. 
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Figure 2. Radio and infrared images of the region surrounding PSR J1622^t950. Each image is centered on the position of PSR J1622^950, indicated by the 
white "+" symbol. Nearby radio sources are labeled A, B, C, and D, and the Hll region, G333.6-0.2, is also indicated. The black contours in panels b and c show 
the MGPS1 0.005, 0.01, 0.03 and 0.05 mJy per beam radio emission levels. The nearby pulsar, PSR J1622^944, is indicated by a black "+" symbol, a) The 
surrounding region as seen in the MGPS 1 radio survey at 843 MHz and at a resolution of 43" . b) The grayscale is the surrounding region as seen in GLIMPSE at 
8 (im at a resolution of 1.2". c) The grayscale as seen by MlPSGAL at 24 (im at a resolution of 6". The white circle indicates the position of IRAS 16190-4946 
with a radius of 23", equivalent to its major axis position uncertainty. 
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Figure 3. Radio and X-ray light-curves of PSR J1622^W50 over 1500 days. The left vertical axis shows the radio flux at 1.4 GHz in units of mJy and the right 
vertical axis shows the absorbed X-ray flux in units of 10~' 3 erg cirr 2 s in the 0.3 - 10.0 keV energy range. The radio flux points are in red where the Parkes 
1 .4 GHz data are denoted by asterisks and the Parkes 6.5 GHz Multibeam pulsar survey measurement is a range of possible values at 1 .4 GHz, calculated from 
a range of possible radio spectral indices, denoted by two crosses connected by a dashed line. (These light-curve points were originally depicted in Figure 1 of 
Levin et al. 2010). The red squares show the ATCA detections of PSR J1622-4950 when extrapolated to 1 .4 GHz. The ATCA detection with the horizontal error 
bar is the average flux value from two ATCA observations, taken in 2009 Dec and 2010 Feb, published by Levin et al. (2010). The ATCA flux errors are the size 
of the data point. The X-ray flux points are in blue where the Chandra observations are denoted by open circles and the XMM observation by a filled circle. The 
error bars are la. The upper-limit on the X-ray flux of PSR J 1622-4950 detected by ASCA in 1999 April (MJD 51291) was F x <3 X 10~ 13 erg cm" 2 s" 1 . 



